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Density  functional  theory  computations  have  been  employed  to  investigate  the  decomposition  of  formic 
acid  (HC02H)  into  C02  and  hydrogen  on  the  |3-Mo2C(101)  surface.  The  adsorption  configurations  and 
energies  of  the  surface  intermediates  (HC02H,  C02,  CO,  H20,  HC02,  C02H,  CHO,  OH,  0  and  H)  have  been 
systematically  characterized.  Among  the  different  dissociation  steps  considered,  our  results  showed  the 
formate  route  (HC02H  ->  H  +  HC02;  HC02  -►  H  +  C02)  is  the  minimum  energy  path  for  hydrogen 
formation  and  C02  has  very  strong  chemisorption.  The  adsorption  and  dissociation  of  formic  acid  on  the 
Mo2C(101)  surface  have  been  compared  with  those  of  Pt  group  metals. 

©  2013  Elsevier  B.V.  All  rights  reserved. 


1.  Introduction 

Selective  and  catalytic  decomposition  of  formic  acid  (FA, 
HCO2H)  into  H2  and  CO2  is  considered  as  one  of  the  potential 
chemical  processes  to  satisfy  the  increasing  energy  demand, 
especially  in  fuel  cell,  green  and  clean  technologies  [1—3].  For  FA 
selective  decomposition  both  heterogeneous  [4—7]  and  homoge¬ 
neous  [8-12]  catalysts  have  been  used.  Recently,  Flaherty  et  al. 
studied  FA  decomposition  on  molybdenum  carbide  and  found  that 
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C-modified  Mo(llO)  and  C-Mo(llO)  are  up  to  15  times  more  se¬ 
lective  than  pure  Mo(llO)  for  H2  formation  [13].  Koos  and  Solymosi 
reported  that  the  highly  stable  molybdenum  carbide  (M02C)  cata¬ 
lyst  prepared  from  the  reaction  of  M0O3  with  a  multiwall  carbon 
nanotube  and  carbon  Norit  can  selectively  decompose  FA  into  CO- 
free  H2  at  373-473  K  [14].  Cui  et  al.  reported  that  Pd  catalysts 
supported  on  M02C  which  is  supported  on  multiwall  carbon 
nanotube  has  much  higher  electrocatalytic  activity  and  stability  for 
FA  electroxidation  than  only  M02C  catalyst  supported  on  multiwall 
carbon  nanotubes  and  only  Pd  catalysts  supported  on  M02C,  and 
they  concluded  that  M02C  is  not  only  a  support  but  also  a  co¬ 
catalyst  [15].  In  addition,  M02C  is  active  for  hydrogenation  and 
dehydrogenation  [16],  as  well  as  for  low  temperature  water-gas 
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shift  (WGS)  reaction  [17,18].  Transition  metal  carbides  which  have 
the  added  benefits  of  lower  cost  with  respect  to  the  Pt  group 
metals,  like  W2C  and  MoC2have  been  found  to  have  Pt-like  reac¬ 
tivity  [19,20]. 

Theoretically,  molybdenum  carbides  have  been  used  to  study 
various  reactions.  The  hydrogenolysis  mechanisms  of  thiophene 
[21]  and  indole  [22,23]  on  clean  P-M02C  have  been  studied  sys¬ 
tematically.  The  chemisorption  and  decomposition  of  small  mole¬ 
cules  [24-29]  such  as  nitrogenous  compounds,  aromatic 
hydrocarbons  and  CO2  were  examined  on  both  0C-M02C  and  P-M02C 
phases.  Theoretical  studies  of  the  chemical  properties  of  methanol 
[30],  methyl  iodide  [31  ],  CO  and  the  promoting  effect  of  potassium 
on  P-M02C  were  reported  systematically  by  Pistonesi  et  al.  [32]. 
Based  on  their  surface  experiments;  and  they  found  that  the 
incorporation  of  potassium  atoms  enhances  the  dissociation  ability 
of  the  C— I  and  C— O  bonds  in  CH3I  and  CH3OH,  while  blocks  the 
dissociation  of  CO.  Tominaga  and  Nagai  built  a  schematic  potential 
energy  surface  for  WGS  reaction  and  concluded  that  C02  formation 
from  CO  oxidation  by  surface  O  is  the  rate-limiting  step  [33].  The 
mechanism  of  CO  hydrogenation  and  the  promoter  effect  of  cobalt 
have  also  been  systematically  reported  recently  [34].  Liu  et  al.  also 
calculated  WGS  mechanism  and  emphasized  the  importance  of 
oxygen  on  the  M02C  surface  [35].  In  order  to  study  the  intrinsic 
WGS  activities  of  M02C,  Schweitzer  et  al.  loaded  Pt  on  M02C  and 
found  M02C  to  play  the  role  of  both  support  and  catalyst  [36].  Shi 
et  al.  [37]  and  Han  et  al.  [38]  calculated  the  surface  energies  of  low 
miller  index  surfaces  of  hexagonal  M02C  to  compare  their  stabilities 
and  concluded  that  the  (Oil )  facet  was  the  most  stable  surface.  The 
elementary  steps  of  syngas  reaction  have  been  systematically 
studied  by  using  ab  initio  thermodynamics  method  by  Andrew  et  al. 
and  M02C  was  proved  to  have  similar  catalytic  properties  with 
noble  metals  [39].  Recently,  Zheng  et  al.  reported  H2  production 
from  ammonia  decomposition  catalysed  by  molybdenum  carbide 
both  experimentally  and  theoretically  [40]. 

In  this  work,  we  carried  out  spin-polarized  periodic  density 
functional  theory  computations  to  study  the  adsorption  and 
dissociation  of  FA  on  the  Mo2C(101)  surface.  Our  goal  is  the  un¬ 
derstanding  into  the  adsorption  configurations  of  FA  and  its 
dissociation  intermediates  on  the  Mo2C(101)  surface  as  well  as  the 
dissociation  paths.  These  results  are  compared  with  those  on  the  Pt 
group  metals,  e.g.;  Pd(lll)  [41-44],  Pt(lll)  [43,45,46],  and  Ir(100) 
[47],  from  recent  computational  studies. 

2.  Computational  details 

2.1.  Model 

M02C  mainly  has  two  crystalline  structures,  the  orthorhombic 
a-Mo2C  phase  [48]  and  the  hexagonal  P-M02C  phase  [49,50].  In  our 
work,  we  used  the  (3-hexagonal  M02C  phase  with  an  eclipsed 
configuration  as  unit  cell  [51,52].  The  calculated  lattice  parameter  of 
the  cell  is  2 a  =  6.075  A,  2b  =  6.069  A  and  c  =  4.722  A,  in  good 
agreement  with  the  experiment:  a  =  b  =  3.011  A  and  c  =  4.771  A 
[53].  Among  all  the  surfaces  of  P-M02C,  the  (101)  surface  was  re¬ 
ported  to  be  most  stable  [37,51-55],  and  there  are  two  types  of  C 
atoms  and  two  types  of  Mo  atoms  on  the  exposed  surface  (Fig.  1). 
For  describing  this  surface  atoms  easily,  the  4-coordinated  (two 
surface  Mo  atoms  and  two  bulky  phase  Mo  atoms)  C  atom  is 
marked  as  CA,  the  5-coordinated  (four  surface  Mo  atoms  and  one 
bulky  phase  Mo  atom)  C  atom  is  denoted  as  Cb.  The  10-coordinated 
(three  surface  Mo  atoms,  three  surface  C  atoms  and  four  bulky 
phase  Mo  atoms)  and  11 -coordinated  (three  surface  Mo  atoms, 
three  surface  C  atoms  and  five  bulky  phase  Mo  atoms)  Mo  atoms 
are  notated  as  MoA  and  MoB,  respectively.  The  total  supercell  con¬ 
tains  a  M032C64  unit  within  a  volume  of  15.40  x  12.11  x  17.96  A,  and 


Fig.  1.  Side  and  top  views  of  the  Mo2C(101)  surface. 


the  exposed  surface  has  16  Mo  atoms  and  16  C  atoms.  In  addition,  a 
smaller  surface  model  containing  4  exposed  Mo  atoms  and  4 
exposed  C  atoms  within  the  red  region  of  the  surface  was  also 
employed  to  study  coverage  effect,  which  is  defined  as  the  exposed 
surface  Mo  atoms,  e.g.;  1/16  ML  for  the  large  model  and  1/4  ML  for 
the  small  model  (Fig.  1). 

2.2.  Method 

All  calculations  were  done  using  the  plane-wave  spin-polarized 
periodic  density  functional  method  (DFT)  in  the  Vienna  ab  initio 
simulation  package  (VASP)  [56-59].  The  electron  ion  interaction 
was  described  with  the  projector  augmented  wave  (PAW)  method 
[60].  The  electron  exchange  and  correlation  energy  was  treated 
within  the  generalized  gradient  approximation  in  the  Perdew- 
Burke-Ernzerhof  formalism  (GGA-PBE)  [61].  For  Mo,  the  core  4p 
states  were  also  taken  into  valence  region  and  totally  12  valence 
electrons  were  included.  The  Kohn-Sham  one-electron  states  were 
expanded  in  a  plane  wave  with  an  energy  cutoff  400  eV,  and  the 
Methfessel— Paxton  scheme  was  used  under  the  electron  smearing 
of  a  =  0.1  eV  [62].  The  vacuum  zone  was  set  up  to  12  A  in  the  z 
direction  to  separate  the  slabs.  The  geometry  optimization  was 
done  when  force  becomes  smaller  than  0.02  eV  A”1  and  the  energy 
difference  was  lower  than  10-5  eV.  For  bulk  optimization,  the  lat¬ 
tice  parameters  were  obtained  by  minimizing  the  total  energy  of 
the  unit  cell  using  a  conjugated-gradient  algorithm  to  relax  the  ions 
and  a  5  x  5  x  5  Monkhorst-Pack  k-point  grid  was  used  for  sam¬ 
pling  the  Brillouin  zone  [63].  The  first  Brillouin  zone  was  sampled 
with  3x3x1  k-point  grid  for  investigating  the  surfaces  (1/16  and 
1/4  ML).  The  nudged  elastic  band  (NEB)  method  was  used  to  locate 
the  minimum  energy  path  [64].  The  computed  vibrational  fre¬ 
quencies  were  used  to  characterize  a  minimum  state  without 
imaginary  frequencies  or  an  authentic  transition  state  with  only 
one  imaginary  frequency. 

The  adsorption  energy  is  defined  as  in  Equation,  Ea dS  =  EA/ 
slab  -  [Esiab  +  Ea];  where  EA/Siab  is  the  total  energy  of  the  slab  with 
adsorbates  A,  Esiab  is  the  total  energy  of  the  bare  slab,  and  EA  is  the 
total  energy  of  free  adsorbate  A  in  gas  phase;  and  the  more  negative 
the  EadS,  the  stronger  the  adsorption.  The  activation  energy  is 
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defined  in  Equation,  Fa  =  Fjs  -  and  the  reaction  energy  is 
defined  as  in  Equation,  Er  =  EFs  -  Fis;  where  E\S,  Efs  and  ETs 
represent  the  total  energy  of  initial,  final  and  transition  states. 

3.  Results  and  discussion 

3.1.  Adsorption  of  the  surface  intermediates 

For  each  adsorbate,  only  adsorption  site  and  geometries  with 
the  largest  adsorption  energy  are  included  while  many  more  have 
been  evaluated.  Table  1  lists  the  adsorption  energies  of  the  most 
stable  adsorption  configurations  on  the  Mo2C(101)  surface  at  1/16 
ML  and  1  /4  ML.  The  rather  small  differences  in  adsorption  energies 
of  these  intermediates  show  that  1/16  ML  is  quite  larger  enough  for 
our  analysis  and  the  coverage  effects  can  be  ruled  out.  Fig.  2  shows 
the  most  stable  adsorption  configurations  of  these  intermediates 
(HC02H,  C02,  CO,  H20,  HC02,  C02H,  CHO,  OH,  0  and  H)  at  1/16  ML; 
and  those  at  1/4  ML  are  given  in  the  Supporting  information.  Our 
discussion  and  comparison  are  mainly  based  on  the  data  at  1/16 
ML,  if  it  is  not  noted  otherwise. 

Two  FA  adsorption  configurations  on  the  Mo2C(101)  surface  are 
located,  the  flat  one  HC02H-flat  (1)  and  the  perpendicular  one 
HC02H— O/OH-down  (2).  As  shown  in  Fig.  2,  HC02H-flat  (1)  has  its 
0=C— 0  group  chelating  with  a  surface  Moa  atom,  and  the  O  atom 
of  the  OH  group  bonding  to  a  neighbouring  surface  MoA  atom,  and 
the  adsorption  energy  is  -1.12  eV,  indicating  a  strong  adsorption.  In 
HC02H-0/OH-down  (2),  FA  stands  over  the  surface  by  the  O  atom 
of  the  0=C  group  interacting  with  one  surface  MoA  atom  and  the  H 
atom  of  the  OH  group  pointing  towards  one  surface  CA  atom;  and 
the  adsorption  energy  is  -1.06  eV,  which  is  only  0.06  eV  less  stable 
than  that  of  HC02H-flat  (1).  This  indicates  that  both  adsorption 
configurations  are  very  close  in  energy  and  they  might  coexist  at 
this  coverage. 

C02  has  two  energetically  very  close  adsorption  configurations, 
C02-0M/CC  (3)  and  C02-0M/0M  (4).  In  C02-0M/CC  (3),  the  O 
atoms  bond  to  surface  MoA  and  Mob  atoms,  and  the  C  atom  bonds 
to  a  surface  CA  atom,  and  the  OCO  angle  is  128.04°;  and  the 
adsorption  energy  is  -0.92  eV.  In  C02-OM/OM  (4),  the  O  atoms 
bond  to  surface  MoA  atoms  and  the  C  atom  bridges  two  MoA  atoms, 
and  the  OCO  angle  is  134.96°;  and  the  adsorption  energy 
is  -0.81  eV.  In  addition,  the  physisorption  is  also  located  and  the 
adsorption  energy  is  only  -0.19  eV,  much  weaker  than  the 
chemisorption. 

CO  has  two  energetically  almost  equivalent  configurations,  CO- 
flat  (5)  and  CO-C-down  (6).  In  CO-flat  (5),  the  O  atom  bonds  to 


Table  1 

The  computed  adsorption  energies  (£ads,  eV)  for  the  most  stable  adsorption  con¬ 
figurations  of  the  intermediates  involved  in  FA  dissociation  at  1/16  and  1/4  ML. 


Species 

1/16  ML 

1/4  ML 

HC02H-flat  (1) 

-1.12 

-1.16 

HC02H— O/OH-down  (2) 

-1.06 

-1.06 

C02-0M ICC  (3) 

-0.92 

-0.87 

C02-0M/0M  (4) 

-0.81 

-0.81 

CO-flat  (5) 

-1.61 

-1.32 

CO-C-down  (6) 

-1.58 

-1.56 

H20  (7) 

-0.66 

-0.63 

HC02— O/O-down  (8) 

-4.04 

-4.10 

HC02-flat  (9) 

-3.44 

-3.54 

HC02— O/H-down  (10) 

-2.88 

-2.76 

C02H-0M/CC(11) 

-2.92 

-2.91 

C02H-0M/CM  (12) 

-2.85 

-2.94 

CHO  (13) 

-3.22 

-3.21 

OH  (14) 

-4.15 

-4.16 

0(15) 

-6.59 

-6.58 

H  (16) 

-2.79 

-2.76 

two  surface  MoA  atoms,  and  the  C  atom  bonds  to  two  surface  CA 
atoms;  and  the  adsorption  energy  is  -1.61  eV.  In  CO-C-down  (6), 
the  C  atom  adsorbs  vertically  atop  on  one  surface  MoA  atom,  and 
the  adsorption  energy  is  -1.58  eV. 

H20  (7)  prefers  the  O  atom  bonding  to  surface  MoA  atom,  and  its 
H  atoms  pointing  towards  a  neighbouring  CA  and  a  neighbouring 
Moa  atom,  and  the  adsorption  energy  is  -0.66  eV. 

The  most  stable  adsorption  configurations  of  formate  are  HC02- 
O/O-down  (8),  HC02-flat  (9)  and  HC02-0/H-down  (10).  In  HC02- 
O/O-down  (8),  HC02  stands  perpendicular  to  the  surface  and  the  O 
atoms  bridge  two  surface  MoA  atoms  with  0-MoA  distances  of 
215  pm;  and  the  adsorption  energy  is  -4.04  eV.  In  HC02-flat  (9), 
HC02  covers  the  surface  and  its  two  O  atoms  bond  directly  to  two 
Moa  atoms  with  the  Mo-0  distances  of  203  pm;  and  the  adsorption 
energy  is  -3.44  eV.  In  HC02-0/H-down  (10),  one  of  its  O  atoms 
bonds  to  a  surface  MoA  atom  with  the  0-MoA  distance  of  207  pm 
and  the  H  atom  points  towards  the  neighbouring  MoA  with  the  H— 
Moa  distance  of  222  pm;  and  the  adsorption  energy  is  -2.88  eV. 
These  show  that  HC02-0/0-down  (8)  is  the  most  stable  adsorption 
configuration. 

Carboxyl  (C02H)  has  two  energetically  very  close  adsorption 
configurations,  C02H-OM/CC  (11)  and  C02H-OM/CM  (12).  In 
C02H— OM/CC  (11 ),  the  O  atom  of  the  0=C  group  interacts  with  one 
surface  MoA  atom  and  the  C  atom  binds  to  one  neighbouring  sur¬ 
face  CA  atom  as  well  as  the  H  atom  points  to  one  surface  CA  atom 
(the  configuration  with  the  H  atom  opposite  to  the  surface  CA  atom 
is  0.03  eV  less  stable,  Supporting  information);  and  the  adsorption 
energy  is  -2.92  eV.  In  C02H-OM/CM  (12),  the  0=C  group  bridges 
two  surface  MoA  atoms  and  the  H  atom  is  opposite  to  the  surface 
(the  configuration  with  the  H  atom  pointing  to  the  surface  is 
0.04  eV  less  stable,  Supporting  information)  and  the  adsorption 
energy  is  -2.85  eV.  These  indicate  their  possible  co-existence. 

Formyl  (CHO,  13)  prefers  the  O  atom  bridging  two  neighbouring 
surface  MoA  atoms  and  the  C  atom  bridging  one  MoA  and  one 
surface  CA;  and  the  adsorption  energy  is  -3.22  eV.  Hydroxyl  (OH, 
14)  has  its  O  atom  bridging  two  neighbouring  surface  MoA  atoms, 
and  the  adsorption  energy  is  -4.15  eV.  Atomic  oxygen  (0, 15)  has 
atop  adsorption  on  one  surface  MoA  atom,  and  the  adsorption  en¬ 
ergy  is  -6.59  eV;  and  atomic  hydrogen  (H,  16)  has  atop  adsorption 
on  one  surface  CA  atom  and  the  adsorption  energy  is  -2.79  eV. 

3.2.  Dissociation  path 

On  the  basis  of  the  above  discussed  adsorption  configurations 
and  energies,  we  analyzed  FA  dissociation  on  the  surface  starting 
from  the  two  most  stable  HC02H-flat  (1)  and  HC02H-0/OH-down 
(2).  The  adsorption  configuration  of  the  transition  states  are  shown 
in  Fig.  3;  and  the  energetic  parameters  are  given  in  Table  2. 

(a)  FA  dissociation:  Since  FA  has  two  adsorption  configurations  in 
close  energy;  we  studied  both  dissociation  routes  accordingly. 
Starting  from  HC02H-flat  (1),  FA  can  dissociate  either  to  HC02  +  H 
(formate)  or  alternatively  to  C02H  +  H  (carboxyl)  or  to  CHO  +  OH 
(formyl  +  hydroxyl).  For  HC02  formation,  the  transition  state 
TS(HC02H-flat/HC02  +  H)  (a)  is  located,  and  the  breaking  O— H 
distance  is  137  pm  and  the  forming  Cb-H  distance  is  also  137  pm. 
The  energy  barrier  is  0.93  eV  and  the  reaction  is  exothermic  by 
1.12  eV.  For  C02H  formation,  the  transition  state  TS(HC02H-flat/ 
H  +  C02H)  (b)  is  located;  and  the  breaking  C-H  distance  is  149  pm 
and  the  forming  CA-H  distance  is  137  pm.  The  energy  barrier  is 
1.47  eV,  which  is  0.35  eV  higher  than  the  adsorption  of  FA 
(-1.12  eV);  and  the  reaction  is  exothermic  by  0.23  eV.  For  FA 
dissociation  into  CHO  +  OH,  the  transition  state  TS(HC02H-flat/ 
CHO  +  OH)  (c)  is  located;  and  the  breaking  C— 0  distance  is  210  pm 
(146  pm  in  HC02H-flat  (1)).  The  energy  barrier  is  0.48  eV  and  the 
reaction  is  exothermic  by  1.16.  It  is  clearly  to  see  that  starting  from 
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276 


HC02H-flat  (1) 
(-1.12  eV) 


219 


CO-flat  (5) 
(-1.61  eV) 


HC02H-0/OH-down  (2) 
(-1.06  eV) 


C02-0M/CC  (3) 
(-0.92  eV) 


HC02-0/0-down  (8) 
(-4.04  eV) 


135 


HC02-flat  (9) 
(-3.44  eV) 


HC02-0/H-down  (10) 
(-2.88  eV) 


C02H-0M/CC  (11) 
(-2.92  eV) 


C02H-0M/CM  (12) 
(-2.85  eV) 


CHO  (13)  OH  (14) 

(-3.22  eV)  (-4.15  eV)  (-6.59  eV) 

Fig.  2.  The  most  stable  adsorption  configurations  of  intermediates  involved  in  FA  decomposition. 


HC02H-flat  (1)  FA  prefers  to  dissociate  into  CHO  +  OH  rather  than 
into  HC02  +  H  or  C02H  +  H  both  thermodynamically  and  kineti- 
cally.  Starting  from  HC02H-0/OH-down  (2),  FA  can  only  dissociate 
into  HC02  +  H  because  of  its  perpendicular  adsorption  configura¬ 
tion.  In  the  transition  state  TS(HC02H-0/OH-down/H  +  HC02)  (d), 
the  breaking  O-H  distance  is  elongated  to  122  pm  and  the  forming 
H-Ca  bond  is  140  pm  (107  and  166  pm  in  HC02H-0/OH-down  (2), 
respectively);  and  the  energy  barrier  is  only  0.02  eV  and  the  reac¬ 
tion  energy  is  exothermic  by  1.18  eV.  This  shows  that  FA  dissociates 
spontaneously  into  HC02  +  H.  Compared  to  the  adsorption  energy 
and  the  dissociation  barriers  of  HC02H-flat  (1),  HC02H-0/0H- 
down  (2)  dissociation  into  HC02  +  H  is  much  more  favourable 
kinetically. 


To  show  this  preference  more  clearly,  we  calculated  the  FA 
adsorption  at  very  high  coverage  (1/2  ML);  and  only  HC02H-0/0H- 
down  (2)  adsorption  configuration  is  possible,  and  the  computed  FA 
adsorption  energy  is  -0.80  eV,  which  is  lower  than  that  of  the  first 
FA  adsorption  at  1/4  ML,  indicating  a  significant  lateral  repulsive 
interaction  of  0.26  eV.  At  such  high  coverage,  HC02H-0/OH-down 
(2)  dissociation  into  HC02  +  H  is  the  only  possible  route. 

(b)  HCO2  dissociation :  As  the  most  favourable  surface  interme¬ 
diate,  formate  dissociation  into  either  O  +  CHO  or  H  +  C02  was 
studied.  For  the  O  +  CHO  route  from  HC02-0/0-down  (8),  the 
transition  state  TS(HC02-0/0-down/0  +  CHO)  (e)  is  located;  the 
breaking  C-0  bond  distance  is  elongated  from  127  pm  to  143  pm 
and  the  energy  barrier  is  2.27  eV  and  the  reaction  is  thermal  neutral 
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TS(HC02-0/0-down/0+CHO)  (e) 
(-1.77  eV) 


TS(HC02-flat/H+C02)  (f) 
(-2.43  eV) 


TS(H+0H/H20)  (i) 
(-0.24  eV) 


TS(C0+0H/C02H-0M/CM)  (j) 
(-2.47  eV) 


TS(C02H-0M/CC/H+C02)  (k) 
(-2.69  eV) 


TS(C02/desorption)  (I) 
(0.00  eV) 


Fig.  3.  The  transition  state  configurations  involved  in  FA  decomposition. 


(-0.02  eV).  For  the  CO2  +  H  route  from  HC02-0/0-down  (8),  it  is 
not  possible  for  a  direct  C-H  bond  dissociation  due  to  its  perpen¬ 
dicular  orientation,  and  the  less  stable  adsorption  configurations, 
HC02-flat  (9)  and  HC02-0/H-down  (10)  from  the  dynamic  bending 
were  considered  [65-69].  Starting  from  HC02-flat  (9)  into  C02  +  H, 


Table  2 

The  computed  energy  barriers  (Ea,  eV)  and  reaction  energies  (En  eV)  for  the 
elementary  steps  of  FA  dissociation. 


Reaction 

Ea/Er 

1/16  ML 

1/4  ML 

HC02H-flat/H  +  HC02  (a) 

0.93/— 1.12 

0.95/ — 1.1 1 

HC02H-flat/H  +  C02H  (b) 

1.47/-0.23 

1.49/— 0.16 

HC02H-flat/CH0  +  CO  (c) 

0.48/— 1.16 

0.49/— 1.12 

HC02H— O/OH-down/H  +  HC02  (d) 

0.02/— 1.18 

0.01/— 1.21 

HC02— O/O-down/O  +  CHO  (e) 

2.27/-0.02 

2.27/-0.04 

HC02-flat/H  +  C02  (f) 

1.01/— 0.26 

1.17/— 0.07 

HC02— O/H-down/H  +  C02  (g) 

0.54/-0.82 

0.54/-0.85 

CHO/H  +  CO  (h) 

1.46/-0.02 

1.53/0.05 

H  +  0H/H20  (i) 

1.36/0.94 

1.31/0.95 

CO  +  0H/C02H-0M/CM  (j) 

1.40/1.02 

1.27/0.89 

C02H-0M/CC/H  +  C02  (k) 

0.23/-0.55 

0.27/— 0.47 

C02/desorption  (1) 

0.92 

0.87 

the  transition  state  TS(HC02-flat/H  +  C02)  (f)  is  located;  where  the 
breaking  C-H  distance  is  146  pm  and  the  forming  Ca-H  distance  is 
174  pm.  The  energy  barrier  is  1.01  eV.  Starting  from  HC02-0/H- 
down  (10)  into  C02  +  H,  the  transition  state  TS(HC02-0/H-down/ 
H  +  C02)  (g)  is  located;  where  the  breaking  C-H  distance  is  147  pm 
and  the  forming  Ca-H  distance  is  162  pm.  The  energy  barrier  is 
0.54  eV. 

Starting  from  HC02— O/O-down  (8),  the  effective  barrier  via  the 
transition  state  TS(HC02-flat/H  +  C02)  (f)  is  1.61  eV;  while  the 
effective  barrier  via  the  transition  state  TS(HC02-0/H-down/ 
H  +  C02)  (g)  is  1.70  eV;  and  the  reaction  is  endothermic  by  0.34  eV. 
Compared  with  the  high  barrier  (2.27  eV)  of  O  +  CHO  formation, 
the  formation  of  H  +  C02  via  the  dynamic  bending  HC02-flat  (9) 
configuration  is  more  favourable  kinetically. 

(c)  CHO  dissociation  as  well  as  H2O  and  CO2  formation:  As  the 
most  favourable  route  for  CHO  +  OH  formation  from  HC02H-flat 
(1),  we  computed  CHO  (13)  dissociation  into  CO  +  H,  and  the 
subsequent  C02H  formation  from  CO  +  OH  as  well  as  H20  forma¬ 
tion  from  H  +  OH.  For  CHO  dissociation  into  CO  +  H,  the  transition 
state  TS(CHO/H  +  CO)  (h)  is  located,  and  the  C-H  distance  is 
189  pm;  and  the  energy  barrier  is  1.46  eV  and  the  reaction  is  almost 
thermal  neutral  (-0.02  eV). 
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Starting  from  the  previously  formed  surface  OH,  the  further 
reactions  of  CO  +  H  into  H20  and  C02H  were  calculated.  For  H20  (7) 
formation  from  H  +  OH,  the  transition  state  TS(H  +  0H/H20)  (i)  is 
located,  and  the  forming  O-H  distance  is  132  pm  and  the  breaking 
surface  C-H  distance  is  135  pm;  and  the  energy  barrier  is  1.36  eV 
and  the  reaction  is  endothermic  by  0.94  eV,  while  the  back  reaction 
has  lower  energy  barrier  of  0.42  eV.  Therefore,  H20  dissociation  is 
more  favourable  than  its  formation. 

For  C02H  formation  from  CO  +  OH,  the  transition  state 
TS(CO  +  0H/C02H-0M/CM)  (j)  is  located,  and  the  forming  C-0 
distance  is  163  pm,  and  the  energy  barrier  is  1.40  eV;  and  the  re¬ 
action  is  endothermic  by  1.02  eV;  while  the  back  reaction  has  lower 
energy  barrier  of  0.38  eV;  these  show  that  C02H  dissociation 
instead  of  formation  is  more  favourable.  Nevertheless,  we  also 
calculated  C02H  dissociation  into  C02  +  H.  It  is  to  note  that  there 
are  two  adsorption  configurations  for  C02H  in  very  close  energy, 
C02H-0M/CC  (11)  and  C02H-0M/CM  (12).  Starting  from  the  more 
stable  C02H— OM/CC  (11),  we  located  the  transition  state  TS(C02H- 
OM/CC/H  +  C02)  (k),  in  which  the  breaking  H-0  distance  is  126  pm 
and  the  forming  H-Ca  distance  is  138  pm;  and  the  energy  barrier  is 
0.23  eV;  and  the  reaction  is  exothermic  by  0.55  eV.  In  addition, 
C02H— OM/CM  (12)  (Supporting  information)  dissociates  into 
C02  +  H  is  less  favourable  kinetically. 

On  the  basis  of  the  computed  energy  barriers  and  the  reaction 
energies  for  C02H  formation  and  dissociation,  C02H  dissociation 
into  C02  +  H  is  kinetically  more  favourable,  while  C02H  dissocia¬ 
tion  back  to  CO  +  OH  is  kinetically  less  favourable,  but  thermody¬ 
namically  more  favourable. 

(d)  CO2  desorption :  Since  C02  has  a  strong  adsorption  on  the 
surface,  the  corresponding  desorption  was  computed.  Starting  from 
the  most  stable  adsorption  configuration  (C02-OM/CC  (3)),  the 
transition  state  TS(C02/desorption)  (1)  is  located;  in  which  the 
distance  of  0-MoA  and  C-CA  is  230  and  227  pm,  respectively;  and 
the  OCO  angle  is  151.36°,  which  is  larger  than  that  (128.04°)  in 
C02— OM/CC  (3);  and  the  desorption  energy  is  endothermic  by 
0.92  eV,  and  this  is  also  the  C02  adsorption  energy. 

The  computed  adsorption  energies  of  the  intermediates 
involved  in  FA  dissociation  in  Table  1  show  negligible  differences 
between  1/16  ML  and  1/4  ML  apart  from  the  flat  adsorbed  CO  (5), 


which  has  a  difference  of  0.29  eV.  Detailed  analysis  into  this  dif¬ 
ference  shows  the  strong  reconstruction  of  the  surface  structure  for 
CO  (5)  at  1  / 4  ML  and  no  such  significant  surface  reconstructions  can 
be  found  for  other  intermediates.  Similar  coverage  effects  have  also 
been  found  for  the  energy  barriers  and  reaction  energies  for  the 
elemental  steps  for  FA  dissociation  (Table  2).  All  these  show  that 
our  surface  models  are  large  enough  for  computing  FA  dissociation. 

3.3.  Potential  energy  surface 

On  the  basis  of  the  computed  energetic  parameters  for 
the  elementary  steps,  we  have  mapped  the  potential  energy  surface 
for  FA  adsorption  and  decomposition  on  the  Mo2C(101)  at  1/16 
ML  (Fig.  4).  FA  has  two  adsorption  configurations  close  in 
energy,  parallel  to  the  surface  (HC02H-flat  (1,  EadS  =  -1.12  eV)) 
and  perpendicular  to  the  surface  (HC02H-0/OH-down  (2, 
fads  =  -1.06  eV)).  However,  they  have  very  different  dissociation 
routes.  Starting  from  the  more  stable  configuration  (1),  a  potential 
route  follows  its  dissociation  (c)  forming  surface  CHO  and  OH 
(CHO  +  OH).  Other  dissociation  routes  towards  surface  formate 
(HC02)  (a)  and  carboxyl  (C02H)  (b)  have  much  higher  barriers  and 
thus  kinetically  much  less  favoured.  Considering  (2),  which  is  a  less 
stable  configuration  than  (1),  the  formation  of  surface  formate  (8) 
proceeds  almost  barrier-less  following  the  dissociation  route  (d). 
We  shall  revisit  this  result  when  discussing  the  selectivity. 

Continuing  the  reaction  pathway  starting  at  (8),  direct  C-H 
bond  cleavage  is  not  possible.  Also  the  dissociation  of  (8)  into  sur¬ 
face  CHO  +  O  is  not  competitive,  as  its  barrier  is  2.27  eV.  Instead, 
dynamic  bending  involving  the  adsorption  configurations  HC02- 
falt  (9)  and  HC02-0/H-down  (10)  become  likely.  These  processes 
need  to  overcome  an  effective  barrier  of  1.70  eV,  and  the  surface 
C02  formation  is  endothermic  by  0.34  eV  (or  1.26  eV  to  gaseous 
C02).  However,  the  overall  reaction  pathway  remains  below  the 
potential  energy  of  the  initial  state  (FA  in  gas  phase  {E  m  0.00  eV)). 

A  potential  competitive  reaction  pathway  may  start  at  the  sur¬ 
face  CHO  (13).  The  dissociation  of  (13)  into  surface  CO  +  H  is 
associated  with  a  barrier  of  1.46  eV,  and  the  reaction  is  nearly 
thermal  neutral.  The  subsequent  elementary  steps  of  CO  and  OH  to 
form  surface  carboxyl  (12)  as  well  as  H  +  OH  to  form  surface  H20 


(CO+H20)  (gas) 


j  0.91 


(8)+(16)  (5)+(16)+(14) 

(13)+(14) 


Fig.  4.  The  potential  energy  surface  for  FA  dissociation  at  1/16  ML  on  the  Mo2C(101)  surface  starting  from  gaseous  FA. 


554 


Q.  Luo  et  al.  /  Journal  of  Power  Sources  246  (2014)  548-555 


Table  3 

Comparison  of  the  adsorption  energies  (£ads,  eV),  dissociation  barriers  ( Ea ,  eV)  and  dissociation  energies  (£r,  eV)  of  FA  on  3-M2C(101),  Pd(lll)  and  Pt(lll)  as  well  as  on  Ir(100) 
surfaces. 


Mo2C(101) 

Pd(lll) 

Pt(lll) 

Ir(100) 

fads  (HC02H) 

-1.12;  -1.06 

-0.39  [41];  -0.62  [42];  -0.40  [43];  -0.40  [44] 

-0.39  [43];  -0.40/-0.34  [45];  -0.26  [46] 

-0.63  [47] 

fads  (HC02) 

-4.04 

-2.37  [41];  -2.71  [42];  -2.52  [43];  -2.49  [44] 

-2.27  [43];  -2.32/-2.50  [45] 

-3.66  [47] 

fads  (C02)a 

-0.92 

0.23  [41] 

fads  (H) 

-2.79 

-2.78  [41];  -2.92  [42];  -2.89  [43] 

-2.78  [43] 

-3.01  [47] 

Ea  (HCO2H  HC02  +  H) 

0.02 

0.58  [41];  1.00  [42];  0.68  [43];  0.49  [44] 

0.72  [46];  0.94/0.88  [45];  0.69  [43] 

fr  (HC02H  HC02  +  H) 

-1.18 

-0.17  [41];  0.00  [42];  -0.27  [43];  -0.09  [44] 

0.03  [43] 

0.18  [47] 

fa  (HC02  -  C02  +  H)b 

1.61 

0.76  [41];  0.88  [42];  1.59  [43];  0.77  [44] 

1.23  [46];  1.56/1.16  [45];  1.88  [43] 

0.87  [47] 

Et  (HC02  ->  C02  +  H) 

0.34 

-0.43  [41];  -0.29  [42];  -0.45  [43];  -0.22  [44] 

-0.43  [43];  -0.31/-0.63  [45] 

-0.06  [47] 

Supercell 

(4  x  4) 

(3  x  3)  [41];  (3  x  3)  [42];  (3  x  3)  [43];  (3  x  3)  [44] 

(2  x  4)  [46];  (3  x  3/2  x  2)  [45];  (3  x  3)  [43] 

(3  x  3)  [47] 

a  Chemisorption  state  energy. 
b  Effective  energy  barrier. 


(7)  need  high  barriers  (1.40  eV  and  1.36  eV,  respectively)  and  are 
endothermic  by  1.02  and  0.94  eV,  respectively.  These  show  that  (12) 
and  (7)  are  metal-stable  intermediates  and  they  prefer  the  disso¬ 
ciation  into  surface  CO  (5)  and  OH.  However,  desorption  of  CO  and 
H2O  is  highly  endothermic  and  would  thus  lead  to  poisoning  of  the 
active  sites.  Since  the  elementary  steps  are  equilibrated,  this  part  of 
the  reaction  may  also  proceed  in  the  reverse  direction. 

With  a  view  to  the  high  selectivity  the  catalyst  provides  for  the 
overall  decomposition  of  FA  towards  CO2  and  H2,  we  should  keep  in 
mind  that  the  entire  potential  energy  surface  remains  below  the 
potential  energy  of  FA;  all  elementary  steps  are  equilibrated. 
Additionally,  we  may  formulate  the  conservation  of  the  number  of 
active  sites  by  the  sum  of  all  intermediates  sticking  to  the  surface, 
which  is  constant.  Due  to  the  thermodynamic  equilibrium,  the 
equation  is  cancelling  down  to  1  =  0Hco2+h  +  ^cho+oh,  where  6  is 
the  coverage.  Although  the  reaction  pathway  via  (13)  also  lead  to 
CO2  formation,  the  probability  is  very  low.  Reason  for  this  obser¬ 
vation  is  the  barrier  of  the  elementary  step  (c),  which  is  higher  and 
requires  a  ~  106  times  longer  lifetime  of  the  transition  state  than 
that  of  (d).  Finally,  the  overall  reaction  rate  for  FA  decomposition  via 
the  HCO2  pathway  was  found  to  be  2300  times  faster  than  that 
following  the  HCO  pathway  (Supporting  information). 

Fig.  4  illustrates  the  potential  energy  surface  and  highlights  the 
most  favourable  route  of  the  decomposition  of  FA  via  formate  to¬ 
wards  CO2  +  2H  (blue  line  in  web  version).  This  agrees  with  the 
experimental  findings  of  CO-free  H2  [14]. 

3.4.  Comparison  with  Pd(lll ),  Pt(lll)  and  Ir(  100) 

Since  molybdenum  carbide  has  been  proposed  to  have  the 
catalytic  activity  of  the  Pt  group  metals,  we  compared  the  energetic 
parameters  of  FA  adsorption  and  dissociation  on  Mo2C(101), 
Pt(lll),  Pd(lll)  and  Ir(100).  Since  the  formate  route  is  the  more 
preferable  mechanism,  it  is  much  easy  for  the  general  comparison, 
and  these  data  are  summarized  in  Table  3. 

As  shown  in  Table  3,  FA  has  much  stronger  adsorption  energy  on 
Mo2C(101)  than  on  Pd(lll),  Pt(lll)  and  Ir(100);  all  these  surfaces 
have  FA  chemisorption.  For  CO2  adsorption,  there  are  quite  large 
differences  among  these  surfaces.  For  example,  CO2  on  Mo2C(101) 
has  much  stronger  chemisorption  energy,  while  Pd(lll)  does  not 
adsorb  CO2  (positive  adsorption  energies). 

For  hydrogen  atom  adsorption,  all  these  surfaces  have  very  close 
adsorption  energies  in  the  range  of  -2.8  to  -3.0  eV;  and  there  are 
no  differences  among  these  surfaces  apart  from  the  fact  that  H 
adsorption  prefers  the  surface  carbon  atom  instead  of  surface  Mo 
atom.  For  formate  adsorption,  M02COOI)  has  stronger  chemisorp¬ 
tion  energy  than  Pd(lll),  Pt(lll)  and  Ir(100). 

For  FA  dissociation  into  surface  formate  and  hydrogen,  Pd(lll) 
and  Pt(lll)  have  very  high  activation  barriers,  while  Mo2C(101)  is 


practically  barrier-less.  Along  with  the  energy  barriers,  FA  dissoci¬ 
ation  into  surface  formate  and  hydrogen  is  strongly  exothermic  on 
Mo2C(101  ),  and  almost  thermal  neutral  on  Pd(lll )  and  Pt(lll ).  It  is 
endothermic  on  Ir(100).  For  formate  dissociation  into  surface  CO2 
and  hydrogen,  Mo2C(101)  has  very  high  effective  barrier,  while 
Pd(lll),  Pt(lll)  and  Ir(100)  have  low  effective  barrier  apart  from 
the  reported  1.59  and  1.88  eV  for  Pd(lll)  and  Pt(lll)  by  Hu  et  al. 
[43],  respectively.  The  corresponding  dissociation  energy  is  endo¬ 
thermic  on  Mo2C(101),  while  strongly  exothermic  on  Pd(lll)  and 
Pt(lll). 

4.  Conclusions 

For  the  catalytic  and  selective  decomposition  of  formic  acid  into 
CO2  and  hydrogen  on  (3-Mo2C(101),  the  adsorption  configurations 
and  energies  of  the  surface  intermediates  (HCO2H,  CO2,  CO,  H2O, 
HCO2,  CO2H,  CHO,  OH,  O  and  H)  have  been  computed  systemati¬ 
cally.  On  the  basis  of  the  most  stable  adsorption  states  the  full 
potential  energy  surface  has  been  mapped  and  the  minimum  en¬ 
ergy  path  has  been  identified. 

Among  the  three  possibilities  for  the  first  step  dissociation  of 
formic  acid  dissociation  (CHO  +  OH,  HCO2  +  H  and  CO2H  +  H),  the 
formate  route  is  most  favourable,  and  surface  formate  from 
HCO2H— O/OH-down  (2)  represents  the  most  important  interme¬ 
diate  and  the  most  stable  adsorption  configuration  HCO2-O/O- 
down  (8)  has  its  two  oxygen  atoms  bridging  two  surface  Moa 
atoms. 

Starting  from  HC02-0/0-down  (8),  formate  dissociation  into 
CO2  and  hydrogen  is  kinetically  much  more  favoured  than  into 
surface  formyl  and  oxygen  (CHO  +  O).  Formate  dissociation  into 
C02  and  hydrogen  is  the  rate-determining  step  and  HC02-0/0- 
down  (8)  represents  the  resting  state.  Our  results  rule  out  the 
possible  formation  of  CO  and  H2O,  and  M02C  exhibits  a  unique 
property  to  catalyse  formic  acid  dissociation  into  CO-free  hydrogen. 

Because  of  the  proposed  Pt-like  properties  for  transition  metal 
carbides,  the  adsorption  properties  of  HCO2H,  HCO2,  CO2  and  H  as 
well  as  their  dissociation  energetic  on  Mo2C(101)  have  been 
compared  with  those  on  Pd(lll),  Pt(lll)  and  Ir(100).  Apart  from 
the  adsorption  of  hydrogen  atom,  which  has  very  close  adsorption 
energies  on  all  these  surfaces,  M02COOI)  can  adsorb  these  in¬ 
termediates  much  stronger  than  Pd(lll )  and  Pt(lll ),  in  particularly 
the  much  stronger  CO2  chemisorption. 

For  the  first  step  of  formic  acid  dissociation  into  surface  formate 
and  hydrogen  on  Mo2C(101),  it  is  barrier-less  and  highly 
exothermic,  while  it  has  higher  barriers  and  is  less  exothermic  on 
Pd(  111 )  and  Pt(  111 ).  Due  to  the  much  stronger  adsorption  of  surface 
formate,  the  dissociation  of  formate  into  surface  CO2  and  hydrogen 
has  higher  effective  barrier,  while  this  dissociation  step  on  Pd  and 
Pt  has  lower  barrier. 
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Appendix  A.  Supplementary  data 

It  is  to  note  that  a  review  article  about  formic  acid  adsorption  and 
decomposition  on  a  set  of  metals  and  metal  oxides  has  been  published 
very  recently.  (http://dx.doi.org/10.1142/S0219633613300012). 

Supplementary  data  related  to  this  article  can  be  found  at  http:// 
dx.doi.org/10.1016/j.jpowsour.2013.07.102. 
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